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ABSTRACT 

Massive evolved stars can produce large amounts of dust, and far-infrared (IR) data are essential for determining the contribution of cold dust 
to the total dust mass. Using Herschel, we search for cold dust in three very dusty massive evolved stars in the Large Magellanic Cloud: R71 
is a Luminous Blue Variable, HD 36402 is a Wolf-Rayet triple system, and IRAS05280-6910 is a red supergiant. We model the spectral energy 
distributions using radiative transfer codes and find that these three stars have mass-loss rates up to 10~ 3 M Q yr _I , suggesting that high-mass stars 
are important contributors to the life-cycle of dust. We found far-IR excesses in two objects, but these excesses appear to be associated with ISM 
and star-forming regions. Cold dust (T < 100 K) may thus not be an important contributor to the dust masses of evolved stars. 
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1. Introduction 

Intermediate-mass asymptotic giant branch (AGB) stars are po- 
tentially the dominant dust source in the Galaxy (Gehrz 1989) 
and in low-metallicity environments like the Large Magellanic 
Cloud (LMC; ISrinivasan et al.ll200 9; Mats uura etai]|2009t) and 
other dwarf galaxies (Boveretal. 2009). However, dust pro- 
duction in high-mass stars (>8 M ) remains uncertain. It has 
been suggested that supernovae (SNe) might be the dominant 
dust factory at high-z, since intermediate-mass stars have not yet 
had time to evolv e into AGB stars (Morgan & Edmunds 2003; 
iDwek et al]|2009l) . However, the amount of dust formed in SNe 
in the local universe is much less than required to expla in the 
dust seen at high-z (e.g., Sugerman et al. 2006; Andrews et al. 
2010). It is also unclear if th e dust fo r ms be fore or after the SN 
explos ion. Alternatively, ISloan et all (J2009) and Valiant e et alJ 
(2009) show that AGB stars can contribute dust at high redshifts. 
These studies point to a need to measure the total dust mass from 
all types of stars to obtain a global picture of dust evolution in 
galaxies. 

Part of the LMC was observed with Herschel (Pilbratt et al. 



2010, this volume) as part of the science demonstration program 
(SDP) and the Legacy program entitled A HERsche l Inventory 
of Th e Agents of Galaxy Evolution (HERITAGE; Meixner et al. 
1201 Ol this volume). In this letter, we describe a first look at 
Herschel Photqdetec tor Array Camera and Spectrometer (PACS; 
Pogli tsch et al.ll2Q10L this v olume) and Spectra l and Photometric 
Imaging REceiver (SPIRE; Griffi n et al.ll2.010l this volume) de- 
tections of 3 examples of dust-producing massive evolved stars 
in the LMC: the Wolf-Rayet (WR) system HD 36402, the lumi- 
nous blue variable (LBV) HDE 269006 (or R71), and the red su- 
pergiant (RSG) OH/IR star IRAS05280-6910. While such stars 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 

** Visiting Scientist at Smithsonian Astrophysical Observatory, 
Harvard-CfA, 60 Garden St., Cambridge, MA, 02138 



have been studied exten s ively in the mid - IR (e.g ., Morri s etaD 
1999 ; IClarketal.1 120031; ICrowtheil 120071: iBonanos et afl T2009: 
Ivan Loon et al.ll2010l) . this study is among the first to probe them 
at A > 200 pm. 

2. Mid- to far-infrared photometry 

We obtained PACS (100 and 160 /an) and SPIRE (250, 350 
and 500 pm) fluxes using apertures roughly the size of the 
source full-width at half-maximum and sky apertures avoid- 
ing regions of high background (Table Q]). We also performed 
aperture photo metry on the Multiband Imaging Photometer for 
Spitzer (MIPS; iRieke et all |2004|) images from the Surveying 
The Agents of Galax y Evolution Spitzer Legacy program (SAGE; 
iMeixner et al1l2006l) . 

In Sect. [3] we examine the optical to far-IR spectral energy 
distributions (SEDs). Optical UBVI and near-IR JHK pho- 
tometry are from the Magellanic Clouds Photometric Survey 
(Zaritskv et al.1 1997D and the Two Micron All-Sky Survey 
(Skrutskie et al1 l2006l) . via the SAGE ca talog. J, K, and L' pho - 
tometry of IRAS05280-6910 are from Ivan Loon et al.l d2005l) . 
/-band photometry for HD 3 6402 is from the D eep Near- 
Infrared Southern Sky Survey (Epchtein et al][l997l) . InfraRed 
An-ay Camera (IRAC; iFazio et all 120041) 3.6 - 8.0 pm pho- 
tometry is from the SAGE catal og. Spectra from th e InfraRed 
Spectrograph (IRS; 5.2 - 38 /urn: iHouck et al.ll2004l) and MIPS 
spect ra (MIPS-SED; 52 -97 pm) are also includ ed (e.g., SAGE- 
Spec; iKemper et al.ll2010i Ivan Loon et alJ 20101). We correct fo r 



extinction using the extinction map from Schle gel et al.l ([1998), 
assuming the stars lie midway the LMC contribution: Ay « 
1 mag for HD 36402 and IRS05280-6910 and A v ~ 0.4 mag 
forR71. 

3. The IR nature of detected sources 

The vast majority of point-sources detected in the LMC SDP 
data appear to be young stellar objects (YSOs) or background 
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Table 1. Target Information and Herschel Flux Densities 



Name 


R.A., Dec. (J2000) 


Type 






F„ (mJy) 












100 nm 


160 yum 


250 yum 


350 yum 


500 yum 


R71 


05 h 02 m 07.39 s ,-71°20'13.1" 


LBV 


447 ±8 


105 + 10 


50 ±6 


49+13 


27 ±9 


IRAS05280-6910 


05 h 27 m 40.11 s , -69°08'04.5" 


RSG, OH/IR 






205 ±8 


46 ± 13 


<11 ± 10 


HD 36402 


05 h 26 m 03.94 s , -67°29'57.9" 


WC4(+0?)+08I 












HD 36402 IR1/YS02 


05 h 26 m 02.99 s , -67°29'57.7" 


YSO? 


531 ±10 


872 + 33 


677 ± 33 


<431±23 


<132 ± 21 



Notes.- Quoted uncertainties are 1 a. Upper limits are either below the 3 a detection level or include substantial flux from an adjacent unresolved 
source. IRAS05280-6910 does not fall within the PACS coverage. The WR star itself is not resolved in Herschel images; it is severely blended 
with an IR source (IR1/YS02). The 350 and 500 pm fluxes of IR1/YS02 include a nearby YSO (N51-YS01), which is unresolved from 
IR1/YS02 at these wavelengths. Aperture corrections for PACS and SPIRE were estimated using the current point-spread functions, and are 
roughly 1.4 and 1.2, respectively. None of the sources are heavily affect by high and/or variable backgrounds (Fig.fT}. 




SPIRE 500 



Fig. 1. Herschel images of R71 (left panels), IRAS05280-6910 
(middle panels), and HD 36402 (right panels, also see Fig. [3]). 
Contours on a linear scale are included where it is difficult to see 
the detection. 



galaxies (Sewi lo et al.ll2010l this volume); AGB stars were not 
detected. Here, we describe the Herschel observations of three 
very dusty massive stars (Fig.[TJ. 

3.1. R71 



R71 (Fig. 0) is an LBV with M birth * 40 M Q dLennon et all 
1993), and is currently experiencing an unprecedented erup- 
tion that began in 2005. One month prior to the SDP observa- 
tions (Oct 2009), it showed a ~2 mag increase in visual bright- 
ness and the V-band light curve was just beginning to plateau 
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Fig. 2. SEDs of R71 (top), IRAS05280-6910 (middle) and 
HD 36402 (bottom), fit to 2-Dust models (orange lines). Far- 
IR excesses are fit to modified blackbodies (green lines). The 
IRS and MIPS-SED spectra are shown in black and the stellar 
components are the dotted lines. In the upper two panels, the 
24 pm points are saturated (lower limits). In the bottom panel, 
the bright far-IR emission originates from IR1/YS02; the WR 
system is detected in PACS, but unresolved from IR1. See text. 



dSzczvgiel et al] l2010ri . Van Loon et al. (2010) find a lack of 
cold dust in R71, as indicated by Spitzer data (T^ < 100 K). 
The Spitzer data were acquired just prior to the current outburst, 
and the Herschel data presented here were obtained at near- 
maximum. However, the PACS and MIPS points appear consis- 
tent with each other, such that a model fit to the optical to MIPS 
data also agrees with the PACS points. This indicates that the in- 
creased emission from the photosphere has not yet significantly 
affected the 100- 160 /mi flux. Note that the MIPS 1 60 pm point 
is an upper limit. 

IRS and ISO spectra of R7 1 show strong PAH, crystalline 
forsterite and enstatite features dVoors et alJll999t iMorris et al] 
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Fig. 3. Composite HST/WFPC2 image of HD 36402 in Ha and 
[S II] (purple). The contours represent (a) IRAC 8 /mi and (b) 
PACS100//m. 



2008> iBuchananet all 120091: 1 Waters! l2010l) . A 10.5 /mi amor- 
phous silic ate feature indicate s domina ntly oxygen-ric h (O-rich) 
chemistry. iMorris et al.l (120081) and lVoors et all (11999) speculate 

that the dust was formed during a prior RSG phase. 

Figure [2] shows R71's SED. Ultraviolet spectra ( Blair et all 

2009) are consistent with the stellar component included in our 
model. We fit the SED with a 2-Dust model (online Table |2j 
lUeta & Meixnerll2003l) . adopting spherical symmetry and grain 
properties of the dominant species: amorphous silicate. Three 
dust components are visible, including a previously unknown ex- 
cess visible in the SPIRE data at 250, 350 and 500 /mi (8, 3.8, 
and 3 cr detections, respectively). The two warmer components 
are modeled by the 2-Dust code using two concentric dust shells. 
We assume th e gas-to-dust ratio for the LMC is if/ - 300 (cf. 
lMeixneretal.ll2010L this volume), yielding M ~ 10~ 6 M yr 1 
for the inner dust shell and M ~ 10~ 3 M yr for the outer shell, 
assuming v^^ =10 kms -1 , consisten t with the wind s peed 
of a shell ejected during an RSG phase. IStahl et all (Il986l) find 
v a 160 km s -1 from the Ha line profile, indicating the MLR for 
the inn er shell could be an o rder of magnitude higher. A dusty 
model dNenkova et al.lfl999l) fits the SED equally well and esti- 
mates a MLR that is the 2-Dust outer shell value. 

We fit the far-IR excess with a modified blackbody: 
F^ oc B^(l d)(l -e~ T| ), where B^(Ja) is the Planck function 
at temperature 7^, r,i is the optical depth, and r^ oc X~P, Here, 
we use B - 1.5. The resulting dust temperature is 9 + 1 K, 
which is extremely low compared to the expected temperature 
of outer shells in evolved stars (~30K; Speck et al. 2000), and is 
instead consistent with temperatures of dense ISM dust clouds 
(see Sect.©. 

3.2. IRAS05280-6910 

IRAS05280-6910 is an RSG OH/IR star. It is detected in all 
Spitzer and Herschel bands up to 350 //m (>3.5 cr; Figs. Q] and 
|2|. The star is heavily extinguished in the optical and near-IR, 
and the 10 and 18 /zm si licate features are seen in absorption in 
the IRS spectrum JKemperet ai]|2010h . Van Loon et al. (2010) 
find no indication of dust colder than ~ 100 K in the MIPS data. 
We fit the SED of IRAS05280-69 10 using 2-Dust (flattened 
geometry) and dusty (spherical symmetry). Both models fit the 
SED reasonably well (only the 2-Dust model is shown in Fig. [2]), 
but predict stronger absorption in the 10 /mi silicate feature 
than is seen in the IRS spectrum, even when including grains 
as large as 1 //m (increasing «o suppresses the silicate feature). 
This discrepancy may be due to the inclusion in the slit of a 
nearby RSG (WOHG347), which is too faint to contribute to 



the SED at other wavelengths, but shows emission near 10 /mi 
with enough flux to veil the silicate absorption by the neces- 
sary amount ( van Loon e t al. 2005). Based on both model fits 
and assuming an outflo w velocity of =s20 k ins' 1 (as measured 
from maser emission; Marshall et al. 2004), the mass-loss rate 
(MLR) was M - 2 - 8 X 10" 3 M Q yr~' when the dust was pro- 
duced (online Tabled. There is no evidence of excess emission 
at A > 100 //m in the SDP data, implying there is no significant 
contribution from cold (<k 100 K) and/or large (/mi - mm size) 
grains. 

3.3. HD 36402 

HP 36402 (Fig. [3]l is a WR star that is part of a triple system 
(Moffat et al. 1990). It is the reddest LMC WR star studied by 
Bonanos et al. ( 2009|) in IRAC, potentia lly due to dust formed in 
colliding stellar winds (ICrowfherll2007l) . HD 36402 is detected 
by PACS (Fig. 3b), which might indicate the presence of some 
cool dust. However, the system is almost totally unresolved from 
the stronger far-IR emission immediately to the west (Figs.[T]and 
|3p), which appears to origina te from a nearby molecular c loud, 
visible in Ha (Fig. [3j Dopit a et al.lll994t IChu et al.ll2005h . We 
refer to this emission as HD 36402 IR1. 

The IRS spectrum of HD 36402 from SAGE-Spec is feature- 
less and reminiscent of R Coronae Borealis stars, which are C- 
rich (no silicate dust) and hydrogen-deficient (no PAHs). The 
continuu m emission is t hought to be due to amorphous carbon 
dust (e.g.. iKraemer et al.ll2005l) . The SED (Fig.© shows a stel- 
lar component and a component from a detached, dusty shell 
(~2 - 24 fim). The far-IR SED is dominated by IR1, so an accu- 
rate fit to the WR system in the far-IR is not possible with this 
dataset. In the lowest resolution images (MIPS 160, SPIRE 350 
and 500 /mi), IR1 is also unresolved from N51-YS01 (upper 
limits in Fig.|2]i, further complicating the SED. A 2-Dust model 
(Fig. |2] orange line) gives M - 7 x 10" 6 M Q yr~'. 

What i s the nature of IR1? Based on its IRAC colors, 
IChu et alj J200l identify it as a YSO (N51-YS02). We have 
attempted to fit a two-component modified blackbody (B - 2.0) 
to the IR1 SED to check if it instead originates from dust heated 
directly by the radiation emanating from HD 36402. These fits 
yield a dust temperature of 64 K, with 15 K dust in its wake, con- 
sistent with this scenario. However, a temperature of 64 K places 
the molecular cloud ~0.3+0.1 pc from the WR system, assuming 
typical luminosities of a WC4 star and two O star companions. 
The high-resolution Ha image (Fig. [3} puts the molecular cloud 
at least 1 .2 pc away from the WR system, assuming a distance of 
50 kpc to the LMC (Schaefer 2008). It thus seems likely that the 
far-IR emission of IR1 is indeed related to a YSO embedded in 
the molecular cloud rather than to direct heating from HD 36402. 
In the HST image (Fig.O, the cloud appears brightest in the re- 
gion adjacent to HD 36402, which may indicate some interaction 
between the two. If this is indeed the case, it is plausible that the 
formation of YS02 was triggered by HD 36402. 

4. Implications 

The new far-IR Herschel data have allowed us to take the first 
steps in assessing the contribution of cold dust to the total dust 
mass in 3 massive stars. For the RSG IRAS05280-6910, the 
models give a total dust mass of *0.3 M . This dust mass is 
extremely large for an RSG star; indeed it is also highly uncer- 
tain, as it is computed from the dust density distribution and the 
grain composition and size distribution, which introduces many 
degeneracies. By adjusting these parameters, we can reduce the 
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dust mass by at least 60%, possibly more. The MLR (online 
Table [2]), while still uncertain, is a more reliable quantity and 
suggests the star has entered a necessarily-brief phase of extreme 
mass loss. The large mass found using our simple, preliminary 
models may indicate a complicated geometry with a preferential 
viewing angle requiring more extensive modeling, or perhaps the 
RSG is embedded in a dusty cloud. Whatever the total envelope 
mass, a lack of excess far-IR emission over the model indicates 
that ancient, cold dust does not contribute significantly. Dust at 
50 K, for instance, must contribute less than 10% of the dust 
mass implied by the model fit to explain this lack of far-IR ex- 
cess. 

The LBV R71 shows three dust components. A component 
emitting at ~10— lOOyum resembles RSG dust. A second compo- 
nent dominates at 3 - 10 //m. Its temperature (T[ n = 490 K) sug- 
gests it formed <50 yrs ago, assuming v ~ 100 k m s" 1 and 7p = 
1000 K, perhaps during the 1970s outburst dWolf & Zickgrafl 
1986). The third, more tentative dust component is visible at 
A > 250 //m. If this feature corresponds to cold circumstel- 
lar dust, the im plied dust mass is Md us t Z 10 _1 M Q , follow- 
ing |Evans_etal. (2003) and assuming the absorption coefficient, 
*-(500yum), is 20 cm 2 g" 1 (lOssenkopf & Hennind[l994T) . This is 
far too much dust for a star of this mass; together with its very 
cold temperature (9 K), this high mass suggests that the far- 
IR emission is pre-existing ISM dust swept up by stellar winds 
and/or is ISM dust along the line-of-sight. Indeed, the contours 
in Figure Q] appear to show diffuse emission at the position of 
R71 at the longest wavelengths. If the far-IR emission instead 
originates from the circumstellar envelope, then very l arge grains 
simila r to the cm-sized grains in the Egg Nebula (IJura et al.l 
2000) might explain the far-IR emission without requiring the 
implied large dust mass and cold temperature. Follow-up spec- 
troscopy or deeper SPIRE imaging may help uncover the nature 
of the far-IR emission. 

The WR system, HD 36402, seems to be forming dust 
in colliding winds, which is too warm to emit much at far- 
IR wavelengths. The apparent far-IR emission from HD 36402 
(Fig. 3b) is unfortunately totally overwhelmed by far-IR emis- 
sion originating from a nearby molecular cloud. Assuming 

7—1 9 -1 

KJlOym) » 140 crrr g and /e(500jum) » 20 cm g 
(lOssenkopf & Hennind Tl994). we find M®jjf = 1.0(+0.3) x 
10~ 2 M Q and M^ - 1.5 ± 0.4 M of dust in the molecular 
cloud. 

Stellar evolution models, while still uncertain, show that 
massive stars like these eventually explode as SNe. Due to its 
proximity to a molecular cloud, the HD 36402 remnant may re- 
semble the SN remnant N49, which has swept up 0.2 M of dust 
from the ISM (Ivan Loonetal.12010 : Otsuka et al. 2010, this vol- 
ume). The absence of strong evidence for very large grains in 
R71 and IRAS05280-6910 does not raise the prospects of RSG 
dust surviving a SN blast. Regardless of their fates, the dust 
masses in these 3 stars are quite large, compared with the dust 
mass found in a typical AGB star, showing that high-mass stars 
are important contributors to the life-cycle of dust even in low- 
metallicity environments like the LMC. However, we emphasize 
that we do not find strong evidence for cold dust and/or large 
grains in any of the three objects discussed here, except where 
it is certain (HD 36402) or likely (R71) to be of interstellar ori- 
gin and not synthesized by the object itself. These observations 
indicate that far-IR data of a much larger sample of luminous 
evolved stars in both Magellanic Clouds will be obtained in the 
full HERITAGE dataset, from which we expect clearer patterns 
to emerge. 
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Table 2. SED Fit Parameters and Results 

T, L, % R, R out /R m R m M total 7^'™ M™^'™ = 
(10 4 K) (1O 5 L ) (K) (R e ) (km) (Af yr~') (K) (Afe) 

~10 u /10" ~10 b /~10 4 9±\ zW 1 

-10" ~10~ 3 

~10 10 7xl0~ 6 
„. „. 64/15 (±1) 1Q- 12 />1.5±0.4 

Notes.- Fitting results from 2-Dust. dusty fits for R71 and IRAS05280-6910 estimate outer shell MLRs 2x to 4x less than th e 2-Dust results. 
IR1/YS02 was fit with two modified blackbodies (64 K and 15 K). F or R71, we use O-deficient silicate optical constants (Ossenkopf et al. 
1992h at long wavelen gths and astronomical silicate optical constants (Draine & Lee 1984) at A < 0.18 jum. For IRAS05280-6910, we use 



R7 1 (inner/outer shell) 


1.5 


4.6 


490/120 


100 


2/1. 


IRAS05280-6910 


0.3 


2.2 


250 


1700 


30 


HD 36402 


18 


4.6 


960 


15 


300 


HD36402IRl(YSO2) 













Ossenkopf et al. ( 1992) O-deficient silicates and for HD 36402, we use amorphous carbon grains from Zubko et al. ( 1996). The 2-Dust models 
use a KMH grain size distribution dKim. MartmT & Hendry 199 4]) with q min - O.Olyum and a = 0.1/mi (a - 1 /mi for IRAS05280-6910). 
Stellar parameters for R71 are consistent with those from lLennon et al] dl993l). r^~ mm is the temperature of the coldest dust component. For 
IRAS5280-6910, u winil is measured from the maser emission (20 kms~'; Marshall et al. 2004). For R 71, we use 10 kms~', but the velocity 
measured from the Ha profile indic ates the inner shel l velocity may be >10x larger (Stahl et al. 1986). For HD36402, a typical velocity for a 
WC4 star is assumed (3000 km s"' : lwillis et alj|2004h . 



